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ABSTRACT 

We examine the damping of non-hnear sub-horizon scale entropy fluctuations in early 
epochs of the universe (T ^ 100 GeV to T ^ 1 keV) by neutrino, baryon, and photon in- 
duced dissipative processes. Results of numerical evolution calculations are presented for 
broad ranges of initial fluctuation amplitudes and length scales. These calculations include 
a detailed treatment of neutrino inflation, neutron and proton diffusion, photon diffusive 
heat transport, and hydrodynamic expansion with photon-electron Thomson drag. Neu- 
trino inflation is treated both in the diffusive heat transport regime where neutrinos are 
optically thick on the length scales of the fluctuations, and in the homogeneous heating 
regime where neutrinos are optically thin on these scales. We flnd considerable conver- 
gence in amplitude evolution for appreciable ranges in initial fluctuation length scales and 
amplitudes. Fluctuations produced with the right characteristics at very early times (T ^ 
100 GeV) are found to survive through the nucleosynthesis epoch. 



Subject headings: cosmology: early universe - hydrodynamics - radiative transfer 
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1. Introduction 

The present model for the history of the early universe envisions significant departures 
from local thermal and chemical equilibrium associated with symmetry breaking events. 
These events might include, for example, the QCD epoch (T pa 100 MeV), the electroweak 
transition (T ICQ GeV), and an inflationary epoch. Though the details and conclusions 
remain uncertain such departures from thermal and chemical equilibrium at early epochs 
have been proposed as production sites for entropy fluctuations. By entropy fluctuations 
we mean any local deviations in the entropy per baryon from the cosmic average value. 
Entropy fluctuations could also result from processes associated with primordial black 
holes, cosmic strings, and domain walls. In this paper we wish to examine how an entropy 
fluctuation produced with a given amplitude and spatial dimension at a given epoch will 
subsequently evolve. We will not be concerned here with how fluctuations are produced, 
but rather with how their amplitudes and spatial scales change in time due to the expansion 
of the universe and neutrino, baryon, and photon dissipative processes. 

The abundances of the elements emerging from the nucleosynthesis epoch can be 
sensitive to the existence of entropy fluctuations (cf. Mathews et al. 1990; Kurki-Suonio 
et al. 1988; Terasawa & Sato 1989abc, 1990; Jedamzik, Fuller, & Mathews 1993; and the 
review of nonstandard Big Bang nucleosynthesis by Malaney & Mathews 1993). Observed 
light element abundances can be used to constrain inhomogeneities at the nucleosynthesis 
epoch. It is hoped that these limits, in turn, may be used to constrain any processes in the 
very early universe which might generate entropy fluctuations. Successfully being able to 
constrain the physics of the early universe in this manner would depend on (1) knowing how 
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fluctuations evolve from production at early epochs through the nucleosynthesis epoch, and 
(2) what the effects of such fluctuations on nucleosynthesis would be. This paper addresses 
the flrst of these issues. 

The production of entropy fluctuations on sub-horizon scales in the non-linear regime 
(amplitude greater than unity) has been discussed extensively in the context of the QCD 
epoch (cf. Witten 1984; Applegate & Hogan 1985; Applegate, Hogan, & Scherrer 1987; 
Kajantie & Kurki-Suonio 1986; Kurki-Suonio 1988; Fuller, Mathews & Alcock 1988; Apple- 
gate 1991; Malaney & Mathews 1993). There is no consensus on fluctuation characteristics 
to be expected from this epoch. In fact, recent lattice QCD results (cf. Brown et al. 1990) 
indicate that the chiral symmetry transition associated with this epoch leads to a second 
order phase transition and, therefore, no phase separation and no fluctuation production. 
However, these results are not definitive, as they are dependent on lattice size and a proper 
treatment of dynamical quarks (Petersson 1991). Fluctuations associated with the QCD 
epoch might be generated by kaon condensate nuggets (Nelson 1990) through a mechanism 
which is independent of the order of the chiral symmetry phase transition. 

The electroweak phase transition is expected to be weakly first order (Kirzhnits Sz 
Linde 1976). It has been proposed recently that baryogenesis may be associated with 
non-equilibrium processes proceeding on the surfaces of bubbles of the low temperature 
phase in the cosmic electroweak transition (Shaposhnikov 1986; 1987; 1988; McLerran 
1989; Cohen, Kaplan, & Nelson 1990; 1991ab; Turok & Zadrozny 1990; 1991; Dine et al. 
1991; McLerran et al. 1991). It may well be that the entropy-per-baryon and, hence, 
the baryon-to-photon ratio resulting from these baryogenesis scenarios, will be left with 
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an inhomogeneous distribution across the horizon (Fuller et al. 1993). The generation of 
these entropy fluctuations is related to the details of the non-equilibrium processes as well 
as the stochastic nature of bubble nucleation and coalescence. 

Other events or processes in the early universe have been suggested as entropy fluctua- 
tion generators. For example superconducting cosmic strings and associated large currents 
and magnetic fields could generate appreciable inhomogeneity in baryon density (Malaney 
& Butler 1989). It is conceivable that an inflationary epoch could generate non-linear 
small-scale isocurvature fluctuations (Dolgov & Silk 1992). For a review of entropy fluctu- 
ation generation mechanism we refer the reader to Malaney & Mathews (1993). 

Our goal in this paper is to describe the evolution in time of the amplitude, size, 
and shape of a non-linear sub-horizon scale fluctuation with properties specified at some 
initial epoch. This problem has been discussed in previous work (Peebles 1965; Misner 
1967; Hogan 1978; Heckler & Hogan 1993; Hogan 1993). Our calculations differ from 
previous studies in that we extend our survey to much higher temperatures (T 100 
GeV), we perform complete numerical multi-spatial-zone calculations including all relevant 
dissipative and hydrodynamic effects, we use a more sophisticated and detailed equation of 
state, and we employ the Boltzmann equation in our calculations of neutrino heat transfer. 
Our results, especially as regards to the survivability of fluctuations produced at very early 
epochs, can differ from previous studies. 

Fluctuations may be produced in the early universe with an adiabatic, isothermal, or 
isocurvature character (cf. Kolb & Turner 1989). It is conceivable that fluctuations could 
be produced which have aspects of all of these characteristics. Adiabatic fluctuations are 
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not fluctuations in entropy-per-baryon, whereas, isothermal or isocurvature fluctuations 
are entropy fluctuations. Once fluctuations enter the horizon they evolve rapidly and 
nearly adiabatically to pressure equilibrium with the background environment (Hogan 
1978). We will term any such fluctuation which has come to pressure equilibrium as an 
isobaric fluctuation. 

In Figure 1 we present a schematic picture for entropy fluctuation classiflcation and 
evolution. In this flgure we show the ratio of entropy density in a fluctuation to average 
entropy density plotted against the ratio of baryon number density in the fluctuation to 
average baryon number density. All fluctuations evolve towards the isobaric line. Hydrody- 
namic evolution tracks to pressure equilibrium for isothermal and isocurvature fluctuations 
are shown as dotted lines. Adiabatic fluctuations cease to exist when pressure equilibrium 
is attained. Once pressure equilibrium for non-adiabatic fluctuations is established subse- 
quent evolution proceeds largely along the isobaric line and is induced by neutrino, photon, 
and baryon dissipative processes. 

In Section (2.) we consider the physics of neutrino, photon, and baryon dissipative 
processes. Section (3.) presents numerical simulations for entropy fluctuation evolution. 
Conclusions are given in Section (4.). 

2. Damping of Isobaric Entropy Fluctuations 

In this section we discuss the evolution of subhorizon scale entropy fluctuations. A 
schematic representation of the kind of fluctuation we consider is shown in Figure 2. In 
this flgure is the mean separation between the centers of fluctuations and L is the 
width of a square wave fluctuation. When fluctuations enter the horizon they undergo 
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rapid hydrodynamic expansion (or contraction) until they come into pressure equilibrium 
with their surroundings. Fluctuations which are initially purely adiabatic are completely 
erased by this process. Fluctuations which have an initially isocurvature or isothermal 
component are not erased by the expansion or contraction to pressure equilibrium. We 
will call fluctuations which have attained pressure equilibrium with their surroundings 
isobaric fluctuations. In this paper we will consider the subsequent damping of isobaric 
fluctuations by neutrino, baryon, and photon dissipative processes. 
We deflne the baryon number overdensity distribution, A(x), by 

ni,{x) ^ nb{l + A{x)), (la) 

where nb{x) is the proper net baryon number density at space coordinate x, and fib is the 
average proper net baryon number density (see Figure 2). We take the average proper 
number density at an epoch with scale factor R to be, 

mix) = m{x){^y , (16) 

where Nb{x) is the average proper number density at T = 100 MeV, and we take R = i?ioo 
when T — 100 MeV. In what follows we choose -Rioo — 1- 

The distribution of energy density in relativistic particles, Sreiix), is deflned in terms 
of the horizon average of this quantity, irei , by 

Erelix) = Srel{l + ^S{x)) . (2) 

In this expression 46 (x) is then the deviation of the energy density from its average value at 
position X. We will be concerned with times early enough that the cosmic energy density 
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is dominated by highly relativistic particles. The average energy density in relativistic 
particles in the dilute high entropy conditions of the early universe is 

erel=9eff{f)af\ (3) 

where a = 7r^/30 , T is the cosmic average temperature, and geff = Ylb9b + i '^fdf 
the number of degrees of freedom in relativistic bosons (6) and fermions (/). In this sum we 
will count only those particles with mean free path, Z, smaller than the spatial dimension 
of the fluctuation. This typical spatial dimension, or fluctuation length, we denote as L 
(see Figure 2). 

In the limit where 5{x) « 1 and in epochs where ^e// c^in be approximated as 
constant in time, we find that 

^, , T{x)-f 

5{x) ^ , (4) 

so that in this limit S{x) is just the deviation from the cosmic average temperature at 
coordinate x. 

We can get a rough idea of how fluctuations come into pressure equilibrium with their 
surroundings if we assume that baryons contribute perfect gas pressure Pt = n^T and 
relativistic particles contribute pressure Prei = \srei- In the following sections we will 
check the validity of these assumptions and we will find significant deviation from perfect 
gas pressure at high temperatures T ^30 MeV. Nevertheless, with these assumptions the 
demand for pressure equilibrium between the interior and exterior of a fluctuation implies, 

^erei{l + ^5{x))^nb{l + l^{x))f{l^5{x)) ^ ^rel + fibf . (5) 
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In the limit where 5 « 1 this expression reduces to 

with s the average entropy-per-baryon in the early universe. Hereafter, a bar over any 
quantity denotes the average value of that quantity over the whole horizon volume unless 
indicated otherwise. 

In the early universe the average entropy-per-baryon is a large quantity, 

s 2.6 X 10^n^^h~^ . (7) 

Here fib is the fraction of the closure density contributed by baryons, and h is the Hubble 
constant in units of 100 km s~^ Mpc~^. The average entropy-per-baryon s is constant 
with time in the limit where we can neglect baryon number violating processes (T ^100 
GeV), black hole production and evaporation, and dissipative processes. Note that s is 
large enough that fluctuations will have 5 < 1 for A < 10^. 

Pressure equilibrium for fluctuations is obtained initially by hydrodynamic expansion. 
The expansion of an entropy fluctuation can be rapid enough that we can regard it as a 
nearly adiabatic process. After pressure equilibrium is established we expect the temper- 
ature in the fluctuation to be lower than T (see Figure 2). Subsequent heat and entropy 
flow into the fluctuation will tend to cause the fluctuation to expand to a new pressure 
equilibrium. Entropy will flow into the fluctuation until the average entropy-per-baryon 
in the fluctuation matches the background average. Heat flow, or other dissipative pro- 
cesses, will of course generate entropy. In practice, however, the entropy generation due 
to dissipative processes is negligible for fluctuations with initial 5 « 1. The fluctuation 



will evolve through a succession of pressure equilibrium states if the heat transport time 
is long compared to a hydrodynamic expansion time. 

The approximate time scale to attain pressure equilibrium between the fluctuation 
and the environment is the sound travel time, Tp, over the fluctuation length L. The 
time scale to transform an initially isobaric fluctuation into an isothermal fluctuation is a 
typical heat transport time scale Th- The ratio of the hydrodynamic expansion time scale 
(or sound crossing time) to the typical heat transport time is, 

Tp L/vs gt I , . 

r2/n ' T ' W 

Th L^/Dh 9effL 

where Vs ~ c/VS is the speed of sound for a relativistic fluid. The heat diffusion constant 
is roughly ~ c{gt/ gef f)li where I is a typical mean free path for heat-transporting 
particles, and gt is the number of degrees of freedom which are effective in heat transport. 
Equation (8a) is only valid in the diffusive limit for heat transport, where I < L. 

Particles with I > L are free streaming on the scale of the fluctuation. Even in this 
limit particles can still transport heat into the fluctuation. For neutrinos with I > L we 
find a heating rate rf^^ ~ {gt/geff)c/l, which yields, 

- ■ ^ (86) 



Qeff I 

We conclude that the assumption that fluctuations evolve through a succession of 
pressure equilibrium conflgurations is valid in the extreme diffusive {I « L) and "homo- 
geneous" {I » L) heating limits. Borrowing terminology from radiation transport studies 
we can describe these limits as "optically thick" and "optically thin" , respectively. 
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For an intermediate mean free path, I r-j L, the heating and pressure equihbrium time 
scales become comparable, and the calculation of the damping of fluctuations by heat 
flow would require a detailed solution of the Boltzmann equation. However, in the case 
of neutrino heat transport in the early universe we expect the approximation of pressure 
equilibrium to hold fairly well, since neutrinos go quickly from optically thick to optically 
thin on the scale of the fluctuation. Furthermore the ratio gt/geff is a small quantity for 
neutrino heat transport in the early universe, which tends to make the ratio in equations 
(Sab) small. 

Self gravity of subhorizon scale baryon number fluctuations is negligible compared to 
pressure forces in the early universe. This is not surprising as the Jeans mass is close to the 
horizon mass in a radiation dominated regime. A rough estimate for the ratio of pressure 
forces Fp to gravitational forces Fg on the edge of a square wave fluctuation gives 



where Lioo is the fluctuation length L measured on a "comoving " scale at an epoch of 
T = 100 MeV in a manner to be described below. 

We will frequently refer proper lengths to a "comoving" scale at T = 100 MeV. The 
fluctuation length, L, and separation length between adjacent fluctuations, L*, as measured 
on the T = 100 MeV scale are denoted by Lioo and Lfpo respectively. At an epoch where 
the temperature is T and the scale factor is R{T) the relation between a proper length L 
and the associated "comoving" length Lioo is 






(10) 
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If the product RT were constant with time then we could write R{T) /R\qq = (lOOMeF) /T. 
However, particle annihilation epochs change the value of RT (e.g. e='=-annihilation; /j,^- 
annihilation; quark-annihilation or QCD-transition) . Figure 3 shows i?(T)(T/100MeV) 
for 10~^ MeV < T < 100 MeV, including a numerical treatment of all relevant annihila- 
tion epochs. In this figure we assume that the quark-hadron transition occurs at higher 
temperatures and we ignore pion degrees of freedom. 

a.) Neutrino Inflation 

The neutrino mean free path Z^, is a rapidly varying function of temperature in the 
early universe, 



where G_p is the Fermi constant, and Ni is the number of relativistic weakly interacting 
leptons at temperature T. These will include , /x^ , , Uf, , u^, , , , u^- , Uj- , 
whenever their masses satisfy mi < T. Similiary Nq is the number of relativistic quark 
flavors at temperature T and may include u , d , s , c , b , and t. For example, at a 
temperature below that for quark annhilation Ni ^ 10 { , /x^, and six neutrino species ) 
and Nq = 0. Here we have neglected neutrino-pion scattering. The dimensionless quantity 
(3 depends on neutrino flavor, the types of relativistic leptons (quarks) in equilibrium at 
temperature T, and the exact lepton (quark) weak couplings. Heckler and Hogan (1993) 
have considered neutrino scattering in a weakly coupled relativistic plasma at zero chemical 
potential. They found that plasma screening corrections are negligible, so that neutrinos 
can be approximated as scattering incoherently off single particles. 




(11) 
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Neutrinos are very efficient heat transporters at any time in the early universe between 
the epochs of the electroweak phase transition (T fti 100 GeV) and neutrino decouphng 
(T 1 MeV). Before the electroweak phase transition neutrino heat transport is inefficient, 
since for vanishing weak gauge boson masses {Mz,w = 0) the neutrino mean free path 
becomes small. Since the neutrino mean free path changes considerably over the history 
of the universe where we follow entropy fluctuations we will have to treat neutrino heat 
transport in two different limits: the homogeneous limit, li, > L; and the diffusive limit, 
< L. 

In the homogeneous limit neutrinos are free streaming through the fluctuation, so 
that neutrinos are optically thin on the scale of the fluctuation L. This implies that 
the neutrino energy density is homogeneous on these scales with e,j{x) = e^- In other 
words, the temperature for neutrinos is now slightly larger than the temperature for other 
relativistic particles inside the fluctuation. In this case neutrino scattering will tend to 
transfer energy into the fluctuation. The heat deposition rate in a fluctuation of size L is 
equal to the number of scattering events per unit time multiplied by the average energy 
transfer per collision. The scattering event rate is the neutrino flux (~ {g,^/lQi)T^) times 
the weak cross section (~ G'pT'^) times the total number of targets within the fluctuation 
(~ (7/8) (2/10) (iV; + 3Ng)T^L^), assuming natural units c=n=l. 

The fractional difference in temperature between the neutrinos and the rest of the 
particles in the fluctuation at position x is just S{x). The average energy transfer per 
collision at position x is then roughly d{x)T^, where Ti, is the homogeneous neutrino 
distribution temperature. Note that T,^ = T. These arguments yield the rate of change of 
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energy density in the fluctuation, 



de{x) 
dt 



FGl5{x)T^ , 



(12) 



where F is a numerical factor. At temperatures T ^ 50 MeV we can approximate F 
{geg^/100) 0.2, where ge = (7/8)2iVe, = 2, and g^ = (7/8)6. This naive estimate for 
the homogeneous neutrino heating rate turns out to be low compared to that derived from 
our detailed numerical treatment. 

We have computed the homogeneous neutrino heating rate by using the Boltzmann 
transport equation. Specifically, we evaluated the collision integral term in the Boltz- 
mann equation. In this calculation we solved for the energy transfer between a relativistic 
neutrino Fermi-Dirac distibution at temperature T and a relativistic Fermi-Dirac distri- 
bution at a lower temperature T{l — 5) . The calculation neglects Fermi-blocking effects on 
scattering particles which is a somewhat reasonable approximation in the early universe. 
Our numerical result has the same form as equation (12) but has F — 1.869, so that the 
net heating rate is nearly an order of magnitude larger than the simple estimate. This 
result is not surprising, however, since the naive estimate of equation (12) does not take 
account of the contribution to heating from annihilation processes such as e'^e~ <-> vP. We 
include these processes in our numerical calculation. These processes dominate the energy 
deposition, accounting for 70% of F. 

In the limit where l^, > L we can employ equations (2) and (3) to rewrite equation 



(12) as 



dt geff lu 



(13a) 
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where, 



- ——F . (136) 



Note that F* is roughly independent of g,^, Ni, and Nq. 

When In < L the optically thick limit obtains and we can approximate the neutrino 
heating as arising from neutrino diffusion. In this limit we can express the time evolution 
of 5{x) for a spherical fluctuation as, 



d5{riQQ) ^ Dh I d f ^ d5{rioo) \ 
dt ~ i?2 rfoo drioo ^ dnoo J 



(14) 



where rioo is the radial coordinate of a spherical fluctuation measured on our comoving 
scale at T = 100 MeV, and D/j ~ {g^/ geff)lv is the heat diffusion constant. For a two-zone 
square wave fluctuation in the limit where l,^ < L we flnd 



| = _A.^i^i^, (16) 
dt geff 



with A* a numerical factor of order unity. In this two-zone approximation for the density 
and temperature distribution of a spherical fluctuation 4S represents the fractional energy 
density difference between the inside (zone 1) and the outside (zone 2) of the fluctuation. 
Matching the heating rates in the diffusive and homogeneous limits for l^, = L we obtain 
A* = F*. 

We now consider a two-zone square wave fluctuation, which will serve to illustrate 
the typical damping histories for fluctuations. A simple prescription for calculating the 
damping of isobaric square wave baryon number fluctuations by neutrino heat transport 
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is as follows. If in time t a fluctuation is heated by dS, it has to adiabatically expand by 
dLiQQ = LiQodS in order to re-establish pressure equilibrium. Thus, 

1 dLioo ^ dS 
Lioo dt ~ dt ' ^ ' 

where Lioo = L/Ras before. Baryon number within the fluctuation is conserved during the 
adiabatic expansion of the fluctuation and so the quantity (1 + A)Lfoo is time independent. 
In order to obtain the rate of increase in radius (or "inflation" rate) of isobaric square 
wave baryon number fluctuations we integrate equation (16), using equation (13) in the 
homogeneous limit, and equation (15) in the diffusive limit. To accomplish this we need 
to know the temperature difference, S, as a function of the baryon number overdensity. A, 
for a fluctuation in pressure equilibrium. 

The pressure due to relativistic particles can be written 

Prel — fp^rel ; (-i'^) 

with fp — 1/3, except during epochs of lepton- or quark- annihilations. Figure 4 shows 
fp{T) during the e+e" annhflation. fp drops from /p = 1/3 to fp = 0.28 at T ^ 0.2 MeV. 
At this temperature a major fraction of kinetic energy of the relativistic gas is converted 
to rest mass energy of e+e~ pairs. This rest mass makes zero contribution to the pressure. 

We consider four temperature regimes in which there are natural limits for the pres- 
sure contributed by baryons. At temperatures below the QCD-transition baryon number 
is carried by color singlet baryonic states. In the regime after the QCD-transition the 
kinematics of baryons can be taken to be approximately nonrelativistic. Below T 30 
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MeV baryons exert perfect gas pressure Pf, 



T ^ nbT{l + A) . 



(18) 



The pressure contribution of the excess electrons needed for charge neutrahty is neghgible 
compared to the baryon pressure. The entropy-per-baryon is so high that the net number 
of electrons is small compared to the number of electron-positron pairs. Requiring pressure 
equilibrium as we did in equation (6), we obtain for T 30 MeV, 



In Figure 4 we show the proportionality constant between A and 6 through the e'^e~ 
annihilation epoch. In order to keep the fluctuation in pressure equilibrium during the 
e+e" annihilation process 5 has to increase. An increase in 5 enhances heat flow into the 
fluctuation. During the e+e" annihilation at T ~ 0.5 — 0.05 MeV there are no particles 
which are efficient heat transporters, so that the pair-annihilation effect is unimportant. 
However, during epochs of quark annihilation at high temperatures or muon and pion 
annhilation at T ~ 100 MeV, the neutrino is an efficient heat transporter and heat trans- 
port efficiency is slightly enhanced during the annihilation epochs. 

Thermally produced baryon-antibaryon pairs are still abundant enough for T ^ 30 
MeV to modify the pressure from the perfect gas limit in equation (18). We find 




(19) 




(20) 



where 



3 




(21) 
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is the total proton (neutron) plus antiproton (antineutron) particle density of a zero chem- 
ical potential Fermi-gas at temperature T. The baryon mass is taken to be m. The 
overpressure exerted by the net baryon number falls below perfect gas pressure, since in- 
creasing the baryon number of a = Fermi gas by two is equivalent to adding one 
baryon and annihilating one antibaryon. For a temperature approaching the Hagedorn 
temperature Th ~ 140 MeV (Hagedorn 1971), the actual baryonic particle density exceeds 
the expression for Upair in equation (21) due to the existence of thermally produced bary- 
onic resonances. Since baryon number can be carried by additional degrees of freedom, the 
resulting overpressure due to extra net baryon number is further decreased for T i— > Th 
(Alcock, Fuller, & Mathews 1987). Note that equation (20) reduces to the perfect gas 
pressure limit of equation (18), when Upair « ni. 

We can estimate a temperature difference 5 between the inside of a fluctuation and 
the cosmic background environment from pressure equilibrium. 



In Figure 5 we illustrate the decrease of fluctuation temperature difference 6 at high tem- 
peratures for three different fluctuation overdensities A. The results in this figure do not 
include baryonic resonances. We show the ratio S/dpg, with Spg the temperature difference 
derived using perfect gas pressure for baryons but with no baryon-antibaryon pairs. It is 
interesting to note that at high temperatures, but after the QCD-transition, fiuctuations 
with net baryon number {ni, ^ 0) can almost coexist in thermodynamic equilibrium with 
a vanishing chemical potential (//^ = 0) phase. 




(22) 
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For times earlier than the QCD-transition baryon number is carried by relativistic 
quarks. The QCD-equation of state in the strong QCD-couphng hmit {as ~ 1) is uncertain. 
The strong couphng limit obtains for T ^ 3-5 GeV (Mueller 1985). In the weak QCD- 
coupling limit (ctg << 1) we can approximate the equation of state as that of a relativistic 
ideal Fermi gas with net baryon number density rib and small baryon chemical potential 
/Xb « T, 

P^ = lN,aT' + 9N-'^ , (23) 

where Ng are the number of relativistic quark flavors, and Pg is the total pressure of quark- 
antiquark pairs plus net baryon number. The overpressure due to net baryon number is 
the second term in equation (23). First order perturbative strong interaction corrections 
would increase this term by a factor (1 — 2as/n)~^ (Mueller 1985). Using equation (23) 
we obtain for the temperature difi^erence between the interior and exterior of a two-zone 
fluctuation 

S = ^—^ — {{A + 1)2 - 1) . (24) 

We are now in the position to integrate equation (16) and flnd by how much a fluctu- 
ation inflates. We take F* — 0.747 and (3 — 2.762 and we assume that these quantities are 
approximately temperature independent. The QCD-transition is taken to be at T = 100 
MeV and the electroweak phase transition is taken to occur at a temperature above T = 100 
GeV. Our estimate treats all quark- and lepton- annihilations as occuring instantaneously 
at temperature T = m. Quark masses are taken here to be = ma = 0, = 100 MeV, 
rric = 1.35 GeV, rrib = 5 GeV and mt > 100 GeV. We take the muon mass as = 105.7 
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and the r-lepton mass as = 1.78 GeV. We take account of the increase of scale fac- 
tor times temperature (RT) and decrease of statistical weight, g, during the annihilation 
epochs. The scale factor R depends on temperature as 

fl=f^yV ^ y\ (25) 

V 10.75/ \lOOMeV J ' ^ ^ 

assuming complete pion and muon annihilation at T = 100 MeV {gg = 10.75). The 
time-temperature relationship is 

for a simple radiation dominated expansion. For the early universe we calculate an average 
net baryon number density, 

fib = 1.7 X 10-^ h?f^. (27) 

In equation (26) we define the total statistical weight in relativistic particles to be = 
Yl:b9b{Tb/TY + {"^ /^)^f 9f{Tf /TY, where we recognize that decoupled particles may 
have temperatures {Tb,Tf) which differ from the plasma temperature (T). The relevant 
statistical weight for entropy density, gs, enters into equations (25) and (27) and is defined 
as gs = j:b9biTb/Tf + (7/8) T.f9f{Tf/Tf. 

We have done numerical calculations for two-zone square wave fiuctuations based on 
the above approximations. We have also generalized this calculation to include many 
zones of different nb{x) and thus obtained a better approximation to the continuum limit. 
We will first discuss our results for two zones and then will present our results for the 
multi-zone calculation in Section 3. 
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In Figure 6 we show the evolution of neutrino mean free path l^P^ = Ij^/R and fluctu- 
ation length Lioo = L/Riov three different overdensities, assuming an initial fluctuation 
length L\qq = 10~^m . The proper fluctuation length of very high amplitude fluctuations 
increases at almost the same rate as the neutrino mean free path until the overdensity A 
is reduced by several orders of magnitude and the neutrinos become optically thin on the 
scale of the fluctuation. Therefore, most of the damping of nonlinear entropy fluctuations 
occurs in the diffusive limit (l^ < L). Note that the neutrino mean free path increase, 
dliy/dt, approaches the speed of light near neutrino decoupling (T fti 1 MeV). The neutrino 
mean free path and the fluctuation length appear slightly discontinous at T = 5 GeV be- 
cause of our oversimplifled treatment of 66-annihilation. However, the discontinuities are 
small because of the large number of degrees of freedom at this epoch in the early universe. 
The results of an exact treatment of annihilation for a single degree of freedom would not 
differ appreciably from our estimates. 

Figure 7 shows the results of our two-zone calculation for the damping of fluctuations 
by neutrino inflation between T = 100 GeV and T = 100 MeV. We display the flnal 
baryon-to-entropy ratio (n^ /s) / of inflated fluctuations as a function of initial fluctuation 
length L\qq for flve different initial baryon-to-entropy ratios {ni,/s)i. We find considerable 
convergence in {nb/s)f for a broad range of initial amplitudes, {nb/s)i, when fluctuations 
have small initial length scales. Fluctuations with very small length scales L\qq 
m will be damped to a characteristic {ni)/s)f ft; 2 x 10~^. Fluctuations with large initial 
length scales show little damping from neutrino inflation. Neutrino inflation becomes an 
important damping process for fluctuations with initial length scales smaller than some 
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critical length. This critical length decreases as the initial fluctuation amplitude decreases. 
In Figure 7 the electroweak epoch horizon scale (T = 100 GeV) is approximately at the 
far right hand end of the horizontal axis. Large overdensity fluctuations with large initial 
length scales have signiflcant damping from neutrinos during the strong QCD-coupling 
epoch. The uncertainties in the equation of state at this epoch translate into uncertainties 
in our damping estimates. We show our results as dashed lines in this uncertain damping 
regime. 

In Figure 8 we display our results for neutrino inflation of fluctuations between T = 100 
MeV and T = 1 MeV. In this calculation we have neglected the effects of baryonic and 
mesonic resonances. Fluctuations with small initial length scales {L\qq ^ 0.1 m) converge 
to a final baryon-to-entropy ratio of {nij/s)f fti 1.1 x 10~^. This corresponds to a baryon- 
to-photon ratio inside the fluctuation of 77 8 x 10~^ during primordial nucleosynthesis. 
Fluctuations with initial length scales larger than {L\qq ^ 1 m) have less damping since 
neutrinos become optically thin on the scale of the fluctuation at low temperatures. The 
dotted line divides fluctuations containing less (left) and more (right) baryonic mass than 
the total baryon mass in the horizon at the QCD-transition epoch (T = 100 MeV). 

The effects of neutrino damping on fluctuations between temperatures of T = 100 GeV 
and T = 1 MeV are shown in Figure 9. This figure combines results from Figure 7 and 8. 
Our calculation assumes no fiuctuation modification from phase transition effects during 
the QCD-transition epoch. The dashed lines indicates uncertainties in the damping due to 
equation of state uncertainties associated with the strong QCD-coupling regime. Even if 
we were to assume no damping of entropy fluctuations in the strong QCD-coupling regime 
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Figure 9 would be essentially unchanged. Significant damping of entropy fluctuations in the 
strong QCD-coupling limit is unlikely because of the small neutrino mean free path during 
this epoch. It is evident that neutrino inflation results in a high degree of convergence 
in the final baryon-to-entropy ratio for a broad range of initial fluctuation amplitudes 
and length scales. If neutrino inflation operated only in the homogeneous limit we would 
expect fluctuations of any initial length scale to converge to the same (nb/s)/. Most of 
the neutrino inflation for fluctuations with length scales in the range 10~^m > L\qq ^10~^ 
m takes place in the homogeneous regime so that we see a fair degree of convergence in 
{nb/s)f for these fluctuations in Figures 8 and 9. 

How sensitive are these results to our assumptions about the microphysical quantities 
P and F*, the expansion rate of the universe, and our simplifled treatment of annihilation 
epochs and phase transitions? When there is sufficient inflation (A/ ^(l/10)Aj with A/, 
Aj the final- and initial- baryon number overdensities) from high initial temperatures (T ^ 
1 GeV) we can obtain an approximate analytic result for the final baryon-to-entropy ratio 

where, 



Quantities in equations (28ab), which are temperature dependent (e.g. g, Nq, and Ni) 
should be taken to be their values at the epoch where the neutrino mean free path is 
of the order of the fluctuation length scale. The (n^/s)/ in equation (28ab) is only very 
weakly dependent on the mean free path (~ /?), the total statistical weight (g), and the 
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number of relativistic quarks Nq and leptons Ni at temperature T. Furthermore, the 
effects of neutrino inflation are clearly insensitive to the average baryon-to-photon ratio 
(~ Vthh'^) and the initial overdensity Aj. There is more significant dependence of the 
final ratio {rn,/s)f on the number of heat transporters {g^) and the efficiency of the heat 
transport (~ F*). Uncertainties in the precise QCD-transition temperature will not affect 
the outcome of neutrino infiation significantly unless the QCD-transition temperature were 
very low (T < 100 MeV). Neutrino inflation is not very efficient after the QCD-transition 
if there are large numbers of thermally produced baryon-antibaryon pairs (Figure 5) as 
would be the case for a QCD-transition temperature T > 100 MeV. We do not expect the 
annihilation of pions and muons to substantially modify our results. 

b.) Photon Infiation and Hydrodynamic Expansion 

At temperatures lower than T 30 keV baryon-to-entropy fluctuations are efficiently 
damped by photon diffusion and hydrodynamic expansion against photon-Thomson-drag 
(Alcock et.al. 1990). The number density of e='=-pairs in the universe remains large com- 
pared to the number density of ionization electrons until the temperature drops below 
about T pa 30 keV. The photon mean free path or, equivalently, the total photon transport 
cross section is determined by photon-electron and photon-positron scattering. When the 
pairs finally annihilate we expect a substantial increase in the photon mean free path. If 
the photon mean free path becomes of order or larger than the fiuctuation size, L, then the 
photons go from optically thick to optically thin on this scale. In the optically thin limit 
all temperature differences between the inside and outside of the fluctuation are erased. In 
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this case we can no longer use radiation pressure to balance the difference in baryon pres- 
sure between the inside and outside of the fluctuation. The unbalanced baryon pressure 
leads to hydrodynamic expansion of the fluctuation. 

In the limit where Klein-Nishima corrections are small (T ^ me) the photon mean 
free path is 

lye ~ ; (29) 

where ut = 6.7 x 10~^^ cm^ is the Thomson cross section, and ngi = rig- + ng+ is the 
number density of electrons plus positrons. At temperatures above the electron rest mass 
the photon-electron scattering cross section is smaller than the Thomson cross section. 
However, the photon mean free path is still small in this case since there are a large 
number of thermally produced charged particles. Figure 10 shows the photon mean free 
path = Ij/R for temperatures below T = 1 MeV. Figure 3 for R{T) can be used to 
convert comoving distances into proper distances. 

In Figure 10 we show the photon mean free path for three different proton fluctuation 
amplitudes Ap. The proton density, Up, which equals the net electron density from charge 
neutrality, is just 

Up = UbAp . (30) 

In Figure 10 we neglect alpha particles. The photon mean free path rises sharply during the 
epoch of e"'"e~-annihilation until the number density of thermal e"''e~-pairs drops below 
the net electron density at temperatures around T 30 keV. After e"'"e~-annihilation 
the photon mean free path depends only on the net electron density and so tracks the 
proton density across a fluctuation. Thus the photon mean free path can vary by several 
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orders of magnitude across a fluctuation. A large photon mean free path contributes to 
efficient damping of smaU-scale nonlinear entropy fluctuations. We consider two limits for 
damping of fluctuations by this process which are the analogues of the two limits considered 
for neutrino inflation. 

Fluctuations expand due to diffusive photon heat transport when the photon mean 
free path is smaller than the fluctuation length scale Ij < L. In this limit the heat diffusion 
equation (14) applies, with the photon heat diffusion constant approximately given by 

~ lyQ 5 (31) 

9eff 

where gt is the statistical weight of heat transporting particles (photons, gt — 2) , and 
gef f is the statistical weight of relativistic particles still "coupled" to the material in the 
fluctuation (7, e^). Using equations (14) and (16) we can obtain an approximate damping 
time scale. We take this to be r^, the time to double the size of a square wave entropy 
fluctuation by photon heat advection, 

1 ^ 1 ' (_L_] ' /^ 2A, + 3A4He \ . . 

^ ~ 5 X lO^s \ m J \lOkeV J \ Ap + 2A^He 7 ' 

where Lioo is the comoving fluctuation length as before, and A4jjg is the '^He density 
deflned in a manner analogous to Ap in equation (30). Since we here deal with temperatures 
well within the nucleosynthesis epoch we must allow for a signiflcant "^He mass fraction. 
Equation (32) assumes complete annihilation of e"'"e~-pairs and thus applies only at low 
temperatures T ^ 20 keV. The denominator in the last factor of equation (32) arises from 
the net electron density within the fluctuation, whereas the numerator is related to the 
temperature deviation, S, in the fluctuation. 
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When the photons become optically thin on the scale of the fluctuation {1^ > L) then 
the temperature across the fluctuation becomes uniform and the fluctuation disintegrates 
by hydrodynamic expansion. The overpressure due to extra baryons and electrons within 
the fluctuation will drive collective fluid expansion out of the fluctuation. However, protons 
moving through an isotropic photon gas will experience a drag force proportional to the 
fluid velocity. As seen from the proton rest frame the photon flirx against the direction of 
motion is slightly higher than the photon flux from the opposite direction. Electrons are 
dragged along by the protons through electric forces. These electrons Compton scatter 
more often in head-on collisions transfering some net momentum to the photon gas. This 
leads to the Thomson drag force on the moving fluid F = {A/3)aTSj{v/c) (Peebles 1971). 
Because of this drag force the expanding fluid will reach a terminal peculiar velocity v = 
(drioo/dt), 

' ' (33) 



where dP /dr iqq is the radial pressure gradient, rioo is the radial coordinate as measured 
on our comoving scale, is the energy density in photons, and rie = n^- — is the net 
electron density. We flnd for the pressure exerted by baryons and electrons below T 30 
keV 

P ^ fn, A» + (n;,,,^ + A,) ' - n;,,, j , (34) 

with the sum running over all nuclei i with nuclear charge Z^. In this expression the 
e"'"e~-pair density divided by fih is 
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which is similiar to equation (21) for baryons. The time scale to double the size of a 
square wave fluctuation by hydrodynamic expansion can be obtained from a characteristic 
fluid velocity and a characteristic fluctuation length Lioo- This time scale is roughly 
Th ~ Liqq/v. After complete e"''e~-annihilation we find 

^ ~ 2.4 X 103s \ m J \10keV J \ Ap + 2A.He J ' ^ ^ 

with the notation as before. This hydrodynamic expansion time scale is similiar to the 
inflation time scale from diffusive photon heat transport in equation (32). It is important 
to note that the damping time scales for diffusive photon heat flow and hydrodynamic 
expansion are independent of fluctuation amplitude. In Figure 11 we show damping rates 
for diffusive photon inflation, (dashed), and hydrodynamic expansion, r^^ (dashed- 
dotted), for a fluctuation with Lioo = 1 m. In this figure we also show the reciprocal of the 
Hubble time r^^ (STrG'e/S)^/^ (solid) between temperatures of T = 100 keV and T = 1 
keV. Damping time scales are shorter than the Hubble time at this epoch for fiuctuations 
smaller than Lioo ^ 1 m and for temperatures below T 30 keV. We conclude that 
small inhomogeneities will disintegrate rapidly below T 30 keV. Larger fiuctuations, 
Lioo > will disintegrate at somewhat lower temperatures. 

In this work we ignore all nonspherical dissipative modes (e.g. convection, dendritic 
instabilities, and Rayleigh-Taylor instabilities). This approximation is valid over all the 
range of conditions we consider except possibly at very late times when rapid hydrodynamic 
expansion against Thomson drag occurs. Simple estimates show that the acceleration at 
these times is so small that the surface of an expanding fiuctuation is Raleigh-Taylor stable. 
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c.) Baryon Diffusion 

High- amplitude subhorizon-scale baryon-to-entropy fluctuations in the early universe 
are modified considerably by baryons diffusing through the primordial plasma (Applegate, 
Hogan, and Scherrer 1987). The baryon diffusion equation in spherical comoving coordi- 
nates, rioo = f/R, is 



where Aj is the ratio of the number density of nucleus i to the total baryon number 
density, and Di is the diffusion constant for the ion corresponding to nucleus i. We include 
in equation (37) a source term representing changes in the nuclear abundances due to 
nuclear reactions, so that the coefficients fij represent nuclear reaction rates. Up to a 
numerical factor of order unity the baryon diffusion constant Di for nucleus i is 



with Vi the thermal baryon velocity, li the baryon mean free path, aij the cross section for 
the scattering of nucleus i on particle A;, and the number density of particle k. 

In Table 1 we list the baryon diffusion constants for neutrons and protons, as well 
as the heat diffusion constant for photons in the early universe (Applegate, Hogan, & 
Scherrer 1987; Banerjee & Chitre 1991; Gould 1993). We show baryon diffusion constants 
for the scattering of baryons on individual particle species k. The effective baryon diffusion 
constant of nucleus i in the primordial plasma due to all scattering is 




(37) 



D, 



1 1 



(38) 





(39) 
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Table 1 also displays the relevant cross sections and the relative cosmological importance 
of individual scattering processes. 

Neutrons diffuse most easily through the primordial plasma. At temperatures above 
T 100 keV neutron diffusion is limited by magnetic moment scattering on electrons and 
positrons. However, at lower temperatures (T < 100 keV), or in high-density fluctuations 
with a local baryon-to-photon ratio ^ 3 x 10~^, neutron diffusion is limited by nuclear 
scattering on protons. Thus, the diffusive damping of high- amplitude fluctuations (A >> 
1) by baryon diffusion is a nonlinear process, since for a higher proton density neutron 
diffusion is less efficient. 

Before the freeze-out from weak equilibrium (T 1 MeV) there is no segregation 
between neutrons and protons, and baryons only diffuse efficiently during the fraction 
of time they spend as neutrons. After weak equilibrium freeze-out protons and neutrons 
diffuse independently, and since protons and neutrons have different diffusion constants the 
result will be spatial segregation of these species. Proton diffusion has negligible effects 
above T ^ SO keV, except for the very smallest fluctuations Lioo << 1"7,. However, such 
small fluctuations are damped efficiently by baryon diffusion before weak freeze-out. 

For temperatures above T fti 40 keV Coulomb scattering of protons off electrons and 
positrons dominates the proton diffusion constant. This is because protons move through 
the plasma as independent particles, in particular protons do not drag along electrons. The 
proton electric charge is shielded by plasma electrons since the Debeye screening length is 
much smaller than a typical interparticle spacing (Applegate, Hogan, and Scherrer 1987). 
The diffusivity of heavy nuclei for temperatures above T 40 keV should as well be 
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limited by Coulomb scattering of the nuclei off electrons and positrons. A simple estimate 
for the diffusion constant of nuclei relative to that of protons can be easily obtained from 
equation (38). The diffusivity of nuclei should be suppressed relative to the diffusivity of 
protons because of a larger Coulomb cross section for the scattering of e+e" off nuclei and 
a smaller thermal velocity of the nuclei. For a nucleus with nuclear charge Zj and nucleon 
number Ai the Coulomb cross section is Zf times larger than that for protons and the 
thermal velocity is 1/Ai times smaller than the proton thermal velocity. This leads to a 
suppression factor of the diffusivity of nuclei relative to that for protons of (l/ZfAi). 

In the temperature range 40 keV ^ T > 0.1 keV Compton scattering of electrons 
off photons limits proton diffusion. In contrast to the situation at higher temperatures, 
the dilute e+e" plasma at lower temperatures is inefficient in shielding isolated charges. 
Charge neutrality then requires a proton and an electron to move together and the effective 
proton diffusion constant is the smaller of the electron diffusion constant and the proton 
diffusion constant. Nuclei have to drag along a cloud of Zi electrons. Naively, we may 
expect the diffusion constant for nuclei to be suppressed by a factor of 1/Zj relative to that 
of the proton-electron system. This is because there are Z^ as many electrons scattering of 
photons and the probability of a change in momentum of the nucleus-electron-cloud system 
is Zi times as large as that for the proton-electron system. Below T 0.1 keV proton 
and nuclei diffusion is again limited by Coulomb scattering as a result of the temperature 
dependent increase of the Coulomb cross section. 

It is instructive to compute the diffusion length of baryons, d{t). This is the root- 
mean-square distance a baryon diffuses in time t. The diffusion length dioo = d/R'va. the 
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early universe becomes (Applegate et al. 1987) 

6 



diooit) = (^-^ D{t')dt'j ' , (40) 

and a simple integration through the e"'"e~-annihilation epoch yields the results displayed 
in Figure 12. In this figure we show dioo neutrons as a function of the temperature of 
the universe for three different local baryon-to-photon ratios (77 = 1 x 10~^, ry = 3.5 x 10~^°, 
and ?7 = 3.5 X lO"-*^^). This figure also shows d^QQ for protons. It is evident that baryon 
diffusion is inhibited in large amplitude fluctuations. 

An approximate time scale to double the size of a square wave fluctuation by baryon 
diffusion is 

- S ■ 

where Dj^ is the relevant baryon diffusion constant and L is the fluctuation length scale. 
Note that the damping time scale is independent of the fluctuation amplitude A. In 
Figure 13 we compare the time scales for neutron diffusion limited by magnetic moment 
scattering (short-dashed), neutron diffusion limited by nuclear scattering in a fluctuation 
with r] = 6 X 10~^ (dashed-dotted) , and proton diffusion (long-dashed) to the Hubble 
time (solid) between temperatures of T = 10 MeV and T = 10 keV. For a fluctuation 
length of Lioo ^ 1 m we flnd proton diffusion to have negligible effects above T pa 20 keV. 
Neutron diffusion is seen to be dependent on the local proton density. A high density 
fluctuation with characteristic interior baryon-to-photon ratio rj ~ 10""^ and fluctuation 
length Lioo ~ 1 m will not be much affected by neutron diffusion whenever the temperature 
is above T 50 keV. 
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3. Evolution of Entropy Fluctuations in the Early Universe 
- a Numerical Simulation 

We have simulated the evolution of high- amplitude baryon-to-entropy fluctuations be- 
tween temperatures of T = 100 MeV and T = 10 keV. The numerical simulation presented 
here does not include the synthesis of elements below temperatures of T pa 100 keV. We 
evolve a spherical region of the universe, which contains an initially Gaussian shape fluc- 
tuation at its center. We assume that the horizon volume is filled with such spherical 
regions so that the mean separation of fluctuation centers is approximatly twice the cell 
radius. These assumptions lead to reflective boundary conditions for baryon number flux 
and heat flux at the spherical cell's edge. We represent a fluctuation by a Lagrangian grid 
with sixty zones in the spherical cell. 

Baryon diffusion, neutrino inflation and photon inflation are treated in a simple in- 
tuitive manner. Initially, we take each zone to be in pressure equilibrium. This pressure 
equilibrium conflguration is characterized by a temperature deviation 5° and proton (neu- 
tron) fluctuation amplitudes (A^) in each zone. Diffusive processes during a time 
interval dt will change the temperature differences in each zone by 6^ i-^ 6^, and the 
baryon densities by A°^^ ^p,n- '^^^^ ^^^^ perturb the existing pressure equality, and 
will result in an adiabatic expansion (contraction) of zones, until pressure equilibrium is 
re-established for new zone values of 6, Ap, and A„. The simple algorithm used here has 
the advantage of being easily tested against our analytic calculations for a two-zone square 
wave fluctuation. 
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To obtain the changes in energy density from diffusive processes in zone i during a 
time interval dt, we rewrite the heat diffusion equation (14) in integral form, which yields 



^-Fy?^i 7 V ~ 3 / X —dx 

R {xi+i-Xi-i) Jo dt 



(42) 



(<-<_i) 

where Xi is the outer boundary of the z-th zone. In this expression Xi is measured on our 
comoving T = 100 MeV scale, so that Xi is short hand notation for x\qq. For homogeneous 
neutrino heat transport we find a change in the energy density in time interval dt given by 

Sl = S^(l-dtF*^^) . (43) 

For this special case the heating is uniform across all zones i in the fluctuation. 

Changes in baryon density due to baryon diffusion are obtained from an integral 
diffusion equation similiar to equation (42). However, efficient neutron diffusion at low 
temperatures necessitates the use of an implicit algorithm. We replaced neutron densities 
at time to on the right-hand side of equation (42) by neutron densities at time ti. The 
resulting matrix equation can be solved by the standard methods of Gaussian elimination 
and back-substitution. In addition we added weak interaction conversion between neutrons 
and protons. Weak reaction rates, as well as the evolution of scale factor R, temperature T, 
and statistical weight g are treated as in the primordial nucleosynthesis code of Wagoner, 
Fowler, & Hoyle (1967) as updated by Kawano (1992). 

Having computed the changes in energy and baryon content of the fluctuation in time 
interval dt, from heat transport, baryon diffusion and weak interactions, we then allow 
each zone to expand adiabatically to a new pressure equilibrium. For this purpose we fix 
the common pressure in all zones after adiabatic expansion to equal the average pressure in 
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the whole spherical cell in which the simulation is performed. Thus, it is easy to compute 
the fractional volume change of zone i, {dV/V)i, needed to attain pressure equilibrium 



We can leave out changes in the baryon density in equation (44) since the associated 
changes in baryon pressure in different zones are negligible compared to the pressure 
changes from heat transport. This approximation is good since the baryon-to-photon 
ratio is everywhere small in our calculation. 

Both photons and neutrinos eventually will decouple on the scale of the fluctuation 
The assumption of pressure equilibrium breaks down when the photons decouple from the 
fluctuation (l^ > L). The simulation of the hydrodynamic expansion of a fluctuation is 
nontrivial since the photon mean free path 1^ varies strongly over the spherical cell radius. 
We encounter situations where photons are still in the optically thick regime in high-density 
regions of the fluctuation but are in the optically thin regime in lower density regions. In 
this case the damping of the fluctuation proceeds by diffusive heat transport in the core 
region of a fluctuation where the baryon density is high, and by hydrodynamic expansion 
near the boundary of the fluctuation. Numerical simulation of this process requires the 
introduction of independent sets of photon zones and baryon zones. 

Neutrinos are easier to treat. The neutrino mean free path li, is independent of the 
baryon overdensity A and neutrinos go quickly from the optically thick limit to the optically 
thin limit. Therefore it is a fairly accurate approximation to switch from diffusive neutrino 
heat transport to homogeneous neutrino heat transport when the neutrino mean free path 
becomes larger than the width a of the Gaussian fluctuation. 




(44) 
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In the following we will summarize the results of our numerical simulation. Figures 14- 
16 show the evolution of three sample high-density spherically condensed Gaussian entropy 
fluctuations for temperatures between T = 100 MeV and T = 10 keV. The initially Gaus- 
sian fluctuations are characterized by four quantities, the average baryon-to-photon ratio 
rj (equivalently Qj^h"^) within the spherical simulation volume, the radius of the spherical 
simulation volume L^qq, the ratio of the radius L|qq to the Gaussian width of the baryon 
density distribution a , and the ratio of the baryon density at the core of the fluctuation 
to the baryon density at the edge of the simulation volume n^, A = {n^ /n^)- For each 
sample we show nine evolution snapshots of the fluctuation at different temperatures. We 
plot the logarithm of the neutron and proton fluctuation amplitudes (Zo^iqA^ and logioAp) 
as functions of the radius rioo = r/R. We reference all proper lengths to the values they 
would have at T = 100 MeV to remove the effects of the overall Hubble expansion. Solid 
lines indicate the local proton fluctuation amplitude, Ap, whereas dotted lines indicate the 
local neutron fluctuation amplitude, A^. 

In Figure 14 the effects of neutrino inflation are evident. Very high- amplitude fluctu- 
ations with an initial central baryon-to-entropy ratio {ni,/s)i ^ 1.5 x 10~^ are damped to a 
characteristic flat-top fluctuation with final baryon-to-entropy ratio {ni,/s)f 1.5 x 10~^ 
(corresponding to an ultimate baryon-to-photon ratio of rjf ~ 10~^). The limiting ampli- 
tude is attained between temperatures of T 100 MeV and T 20 MeV depending on 
initial fiuctuation parameters. 

The modifications of fiuctuations due to neutron diffusion are evident in all three 
sample evolutionary plots. Figure 16 illustrates that small-scale fiuctuations, L ^ 0.5 
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m, are almost completely homogenized and damped out before the onset of primordial 
nucleosynthesis at T a; 100 keV. 

Photon inflation and the hydrodynamic expansion of fluctuations at low temperatures 
are seen to be rapid and rather violent events. This is evident in Figure 15. Between 
temperatures of T 35 keV and T pa 10 keV these processes lower the core proton density 
of fluctuations by several orders of magnitude. The effects of proton diffusion at low 
temperatures (T fti 20 keV) add to the violent decay of such small-scale fluctuations. 

4. Conclusions 

In the present study we have investigated the evolution of fluctuations in the baryon- 
to-photon ratio (entropy fluctuations) between the end of a cosmic electroweak phase 
transition at T f» 100 GeV and the end of primordial nucleosynthesis at T 5 keV. Our 
study focused on nonlinear sub-horizon scale fluctuations. Entropy fluctuations may result 
from an epoch of the early universe where there is a departure from local thermodynamic 
equilibrium such as a flrst order electroweak phase transition. We have shown that entropy 
fluctuations evolve rapidly to an isobaric character. 

Fluctuations are found to damp by five different physical processes: neutrino inflation 
in the diffusive and homogeneous limits; photon inflation and hydrodynamic expansion; 
and bar yon diffusion. 

Neutrino inflation provides the dominant damping of entropy fluctuations between 
temperatures of T fa 100 GeV and T 1 MeV. Fluctuations are damped by neutrino 
heat transport in such a way that almost any initial fluctuation amplitude converges to 
a generic flnal fluctuation amplitude. Fluctuations of high initial amplitude and with an 
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initial fluctuation length Lioo ^ 10~ m become damped to a characteristic baryon-to- 
entropy ratio {ni,/s) fti 2 x 10~^. These fluctuations evolve to this amplitude between 
T Ri 100 GeV and T fa 1 GeV, depending on their initial characteristics. Fluctuations in 
the length scale regime 10~^ m ^ Lioo ^ 10~^ m are damped to a characteristic baryon- 
to-entropy ratio (n^/s) ~ 1.1 x 10~^ (corresponding to an ultimate baryon-to- photon ratio 
r]^8x 10-^) at an epoch between T 100 MeV and T 1 MeV. 

The effects of photon inflation and hydrodynamic expansion on entropy fluctuations 
are shown to be significant at temperatures below T 30 keV, corresponding to the ap- 
proximate completion of e"'"e~-annhilation. Fluctuations with an initial fluctuation length 
Lioo ^ 1 m are almost completely erased by these processes by the time the temperature 
has reached T f» 20 keV. 

Our numerical calculations of fluctuation evolution include all relevant diffusive and 
hydrodynamic damping processes. We do not include any modifications due to nuclear 
reactions here. Clearly our results have important implications for primordial nucleosyn- 
thesis yields in inhomogeneous cosmologies. We will address these issues in future work. 
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Figure Captions 

Figure 1: We show a schematic representation of the ratio of entropy density in fluctua- 
tions to the average entropy density plotted against the ratio of the baryon number den- 
sity in fluctuations to the average baryon number density. We show adiabatic, isothermal, 
isocurvature, and isobaric (heavy line) fluctuations. The dotted lines represent the tracks 
for the nearly adiabatic hydrodynamic expansion/contraction of fluctuations to pressure 
equilibrium. 

Figure 2: Three square wave fluctuations are shown. The net baryon number density 
ni,{x) and the plasma temperature T{x) is shown as a function of length scale x. The size 
of the high density region of a fluctuation, or fluctuation length scale, is L, and the mean 
separation between centers of fluctuations is L*. 

Figure 3: The product of cosmic scale factor and temperature R{T){T /lOOMeV) plotted 
against temperature T in MeV. We take i? = 1 at T = 100 MeV. 

Figure 4: The proportionality constant fp (solid line) between pressure exerted by rela- 
tivistic particles and the energy density in relativistic particles plotted against temperature 
T in MeV. Also shown is the temperature dependence of the proportionality constant be- 
tween the deviation in relativistic energy density and the baryon overdensity within a 
fluctuation (dashed line). 

Figure 5: We plot the ratio of the temperature (or energy density) deviations in a fluc- 
tuation for two different baryon equation of state against temperature. Here Spg is the 
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deviation using a perfect gas Maxwell-Boltzmann equation of state for neutrons and pro- 
tons, while 6 is the deviation including baryon-antibaryon pairs. We show this ratio for 
three baryon overdensities A. 

Figure 6: The proper radius in units of 10~^ m multiplied by is shown as a function 
of temperature for three different initial fluctuation amplitudes. In all three fluctuations 
an initial radius Lioo = 10~^ m is assumed. The fluctuations are distinguished by different 
initial baryon-to-entropy ratios {ni,/s)i = 10~^, 10"'^, and 10~^. The neutrino mean free 
path is also displayed (dashed line). 

Figure 7: Displayed are the effects of neutrino inflation between the epochs of T = 100 
GeV and T = 100 MeV on entropy fluctuations. We give the flnal baryon-to-entropy ratio 
(n^/s)/ of neutrino-inflated two-zone fluctuations as a function of the initial fluctuation 
radius L\qq in meters. We show five different initial (T = 100 GeV) values of {ni)/s)i. 
Dashed lines indicate that {ni,/s)f depends on uncertainties in the QCD equation of state 
in the strong coupling limit. The approximate horizon scale at T = 100 GeV (denoted 
electroweak horizon EWH) is on the far right hand side of the horizontal scale. 

Figure 8: Effects of neutrino inflation between the epochs T = 100 MeV and T — 1 
MeV. Notation is as in Figure 7. We give (rib/s) / for fluctuations with initial length scale 
Lj^oo in meters and initial (T = 100 MeV) amplitude {nb/s)i as labeled. The approximate 
horizon scale at T = 100 MeV (QCDH) is shown. The dashed line divides fluctuations 
containing more (right side) and less (left side) baryonic mass than the baryon mass in the 
QCD-horizon. 
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Figure 9: Effects of neutrino inflation on entropy fluctuations over tlie liistory of the 
universe between epochs of T = 100 GeV and T = 1 MeV. Notation is as in Figures 7 and 
8. 

Figure 10: Photon mean free path = l^/R in meters as a function of temperature. 
Results are displayed for three different proton fluctuation amplitudes Ap. A fraction of 
the closure density in baryons Qb = 0.013/j.~^ has been assumed. 

Figure 11: The fluctuation damping rates due to diffusive photon inflation (dashed) and 
hydrodynamic expansion against Thomson drag (dashed-dotted) and the Hubble expansion 
rate (solid) are given as functions of temperature. The fluctuation length is taken as 
Lioo = 1 m. 

Figure 12: The proton- (dashed) and neutron- (solid) diffusion lengths dioo = d/R in 
meters given as functions of temperature. We assume that the diffusive random walk starts 
at T = 100 MeV. The neutron diffusion length is shown for three different baryon-to-photon 
ratios rj. 

Figure 13: The inverse of fluctuation damping time scales from baryon diffusion are 

displayed as functions of temperature. Shown are the neutron diffusion damping rate 

limited by ne-scattering (short-dashed), neutron diffusion damping rate limited by np- 

scattering (dashed-dotted), and proton diffusion damping rate limited by e7-scattering 

(long-dashed). For comparision the Hubble expansion rate (solid) is given. The key in the 

upper left-hand corner of the flgure displays the dependence of the different time scales 
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on assumed fluctuation radius Lioo = 1 m, average fraction of closure density in baryons 
rib = 0.25h~^, and proton fluctuation ampUtude Ap = 1. 

Figure 14: The evolution of a sample entropy fluctuation over the history of the universe 
between epochs T = 100 MeV and T = 30 keV. We show nine evolution snapshots at 
different temperatures T. Each evolution snapshot displays the logarithm of the proton 
fluctuation amplitude logioAp (solid) and the logarithm of the neutron fluctuation am- 
plitude Zo^ioA„ (dotted) as functions of radial coordinate rioo = r/R in meters. The 
assumed initial fluctuation characteristics are displayed at the bottom of the plots. 

Figure 15: Same as Figure 14 but with different initial fluctuation characteristics as 
shown. 

Figure 16: Same as Figure 14 but with different initial fluctuation characteristics as 
shown. 
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